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Abstract 

We numerically study a triangulated surface model in by taking into account 
a viewpoint of string model. The models are defined by a mapping X from a two- 
dimensional surface M to , where the mapping X and the metric g oi M are the 
dynamical variables. The sum over g in the partition function is simulated by the 
sum over bond lengths and deficit angles by using the Regge calculus technique, and 
the sum over g is defined to be performed independently of the sum over X. We 
find that the model undergoes a first-order transition of surface fluctuations, which 
accompanies a collapsing transition, and that the transitions are reflected in the 
internal geometry of surface. Fluid surface models are also studied on dynamically 
triangulated surfaces, and the transitions are found to be of second order. The order 
of the transition remains unchanged from that of the conventional model defined 
only by the variable X both in the fixed-connectivity and the fluid models. 
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1 Introduction 



Two-dimensional surfaces are interesting objects in the sense that the length 
scales are ranging from microscopic scales to macroscopic ones. The micro- 
scopic string pp and the macroscopic biological membranes P|5f?|f5] can be 
described and hence unified by a surface model of Helfrich and Polyakov [BIT]. 

The models are defined by a curvature Hamiltonian, which is given by an ex- 
trinsic curvature energy p[7|8] . The surface shape or its motion is considered 
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to be governed by the curvature Hamiltonian. One interesting phenomenon 
is a phase transition of shape transformation, which separates a flat phase 
at large bending rigidity from a crumpled phase at small bending rigidity 
PllUfllfl2fl3fl4] . Numerical studies have been devoted to understand the 
phase transitions [T5|16|17j . While recent simulations show that the model 
undergoes a first-order transition [T8|19|20|21|22j . a continuous transition is 
also predicted by theoretical studies based on the renormalization group tech- 
nique (HHl]. Thus, the transition still remains to be studied. 

In the string model context, a two-dimensional surface M and the image X{M) 
of the mapping X from M to R"^ are independently treated. The length scale of 
M is not always identified to the one in X{M). This is because the mapping X 
and the metric g of M are considered as the two different dynamical variables. 
To the contrary, in the case of numerical studies for real membranes, the two- 
dimensional surface M is always fixed and hence is not taken into account to 
define the model. This corresponds to the case where the Euclidean fiat metric 
gab = Sab is assumed in M. In this case, the length scale in M is fixed; the bond 
length of the triangulated M is fixed to some constant. In the case of the 
induced metric gab = daXf^dbX^, M can be identified with X{M). In this case, 
the length scale of M is identified to that of the external space R^, which is 
also an Euclidean fiat space. The surface position X{M) is the only dynamical 
variable in both cases. Thus, little is known about the dependence of the phase 
structure on the variable g in the triangulated surface models. Therefore, it 
is still interesting to study numerically the surface model described by both 
of the variables X and g and to see whether the phase structure of the model 
depends on g or not. 

There have been a lot of studies on models defined by a mapping from a 
D-dimensional surface M to the d- dimensional space R'^, including d = 
the matrix model [23], from the viewpoint of string model [21], where the 
mapping X, including M, is considered to be dynamically generated. However, 
the phase structure of surface models has not been so extensively studied at 
least numerically by assuming the metric as a dynamical variable in the low- 
dimensional cases d = 2,3. 

In this paper, we numerically study a surface model defined by a mapping X 
from M to R^, where the simulation is computationally less time-consuming 
than R^. The variable g is summed over by using the Regge calculus technique 
[25|26|27j in the partition function. M is topologically fixed as a sphere, and 
as a consequence the sum over topology is not included in the sum over g. 

This paper is organized as follows: in Section [2] the model is defined on both 
fixed-connectivity surfaces and dynamically triangulated surfaces. In Section 
[3] we describe how to obtain the model from the continuous Hamiltonian for 
strings. A discrete metric tensor is introduced to obtain the discrete model. In 
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Section H] detailed information of Monte Carlo (MC) technique is given. The 
numerical results including those of the conventional models are presented in 
Section O and we summarize the results in the final section El 



2 Model 



{D = 2) 



M 



{d=2) 




(a) 




Fig. 1. (a) The tangential plane ^{^) of (b) a triangle A in M, and (c) the image 
X{A) in R'^ (d=2), where ^' is a coordinate mapping from A to R^(Z) = 2), and X 
is a mapping from M to R'^ {d = 2). The triangle Z\ in M is almost flat but slightly 
curved, while ^{A) and X{A) are assumed to be linear. 

Let M be a two-dimensional spherical surface, which is not always included 
in the three-dimensional Euclidean space R^, and M is assumed to be trian- 
gulated with smooth triangles. The reason why M is triangulated is because 
the numerical studies including the one in this paper are always performed on 
triangulated surfaces. The triangulation is characterized by the three numbers 
A^, A'^b(= 3A^ — 6), and Nt{=2N — 4), which are the total number of vertices, 
the total number of bonds, and the total number of triangles, respectively. A 
linear triangle in Ti^{D = 2) corresponds to a triangle in M by a coordinate 
mapping \E'; the triangle edges are considered as local coordinate axes on M. 
M is also assumed to be embedded in the external space R'^ ('^ = 2) as X{M) 
by a mapping X. Fig|T]shows the linear triangle "^{A) in R^(D = 2), the trian- 
gle A in M, and the image X{A), which is also a linear triangle in R'^ {d = 2). 
If M is included in R'^((i = 2), the triangle A in Fig{Tt^b) and the tangential 
triangle "^{A) in Fig. [It^a) should be identified with the triangle X{A) in Fig. 



Our basic assumption is that M is not always included in R'^((i = 2). In order 
to define a length scale in M, we identify the edge length L of Z\ with that 
of "^{A). It is not unreasonable to identify L of Z\ with that of "^{A) because 
the length scale in M can locally be fixed by a coordinate mapping \E'; where 
the triangle surfaces are not always identified between A and "if {A). Thus, the 



He). 



3 



triangle inequalities 

Li + Lj > Lk (1) 

are satisfied on A in M. Only difference between A and '^{A) is in the internal 
angles and in the area of Z\, where is defined by 

^A = 7: (^1^2! sin $3 1 + L2L3I sin$i| + L3L1I sin $2!) , (2) 


where |sin$j| is because $j is not always constrained to be 0<$j<7r. While 
the sum of internal angles of the linear triangles is given by Y^\=i = ^ 
(0 < $° < tt), and Yh=i 0i = (0 < 0i < tt), the sum Yl\=i is not always 
identical to vr on Zi. The condition < $j < vr corresponds to the ordinary 
triangles such as the one shown in Fig. EJ^b), and moreover $j [O?^]) ^Iso 
corresponds to those triangles since | sin$j| in Eq.([2]) is defined to represent 
the area of those triangles. 

The deficit angle 99 of Z\ is defined by 

^ = j2^,-n. (3) 

1=1 



Therefore, the internal angle of A can be expressed by 



$, = $ni + -), = 1,2,3). (4) 



TT 



We should note that $° is given only by the edge length L on '^{A). Since 
L is identified to the one of A, $j is given by using only L and if. Thus, the 
sum over L and ip can simulate the sum over metric g on M, where and L 
are assumed to be independent of each other on A. This is a Regge calculus 
approach to the sum over metrics g on M [27j . 

It is possible that \\mj^.^{l / Nt) Yh=i V5j = is violated: 

Ito^Ev'./O. (5) 

The reason of this is as follows: Y^f=i (fii is not always identical to the sum over 
deficit angles Y^f=i Si, where 5i is the deficit angle defined by 5i = 27r—J2j{i) 
where is an internal angle of the triangle meeting at the vertex i. 
In fact, if M is piece-wise linearly triangulated, we have J2iLi = 2vrx = An 
while Y^^Ji (fi = 0, where x is the Euler number. On the contrary, we assume 
that M is smoothly triangulated, where "smoothly triangulated" means that 
every Z\ is a smooth triangle and 5j = at every vertex i. In this case we have 
Yl^iLi Si = while ipi = 47r, which is the prediction of the Gauss-Bonnet 
theorem / ^d?xK = 2ttx = 4vr on smoothly triangulated M, where K is the 
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Gaussian curvature. Therefore (1/A^t) Y^iJi V^j = (N^oo) is satisfied if M is 
piece-wise linearly triangulated or smoothly triangulated. However, the above 
mentioned basic assumption for M does not always imply that M is either 
hnearly triangulated or smoothly triangulated, because the constraints (pi = 
and 6i = are not imposed on the triangulation. 

We comment on the range of (p. A constraint on R) is given by the 
integration measure described below in this section. As a consequence, Iv^l is 
limited to have a value in some finite range in R. Thus, the internal angles 
of a triangle i are automatically determined by Eq. (jlj) from a given (pi. 
Therefore $ is not always constrained to be < $ < vr as mentioned above. 
Informations on the value of ip, the variance of Z^^^i (pi, and the range of ip 
are given in the final part of Section [5l 

Although the coordinate axes are assumed to be varied independently of the 
mapping X, the dynamical triangulation technique is still interesting from the 
viewpoint of reparametrization invariance in the discrete model. The model 
on dynamically triangulated fluid surfaces is also studied in this paper. 

The discrete Hamiltonian S of the model on the triangulated surface M is 
defined such that 



S{X,{L,^}) = S, + bS2, 

^i = ^E^i(^)M^' S2 = ^Y.S2iA)/A^, (6) 

A ^'^ A 

where 5*1 is the Gaussian bond potential, and 5*2 is the bending energy. in 
5*1 and 5*2 is the area of the triangle Z\ in M and is defined by Eq.([2]). S\{A) 
in 5*1 is defined by 



5i (Z\) = i\ [Li + Li) - 2iii2 cos 03 L1L2 cos $3 
+£l [lI + Lj) - cos 01 L2L3 cos $1 

+il (lI + Ll) - 2£s£i cos 02 L3L1 cos $2, (7) 



where the symbols Lj, $j are the bond length and the internal angles of A 
in M, and £j, 0j are those of X{A) in R°'((i = 2). The symbol ^2 {A) in ^2 is 
given by 



^2 (A) = 2L\ (1-no-ni) - LiL2Cos$3 (ni -no) ■ (n2-no) 
+2L2 (l-no-n2) - L2L3Cos<l>i (n2-no)-(n3-no) 
+2L3 (l-no-n3) - L3LiCOS$2 (n3-no) ■ (rii -no) 
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where no (g Z2 = {— 1, 1}) is a unit normal vector of the triangle X{A), and 
n, is a unit normal vector of the nearest neighbor triangle i. The value of n 
is naturally defined by the orientation of the triangle. Figure |2]|^a) shows two 
different values of n, and Fig. [2](b) shows Hq and rij (i = l,2,3) corresponding 
to those in S2 (Zi) of Eq.®. We should note that n ^ R^, and hence the model 
is defined within R'^((i = 2). 



C R2 





(a) (b) 

Fig. 2. (a) Two possible values of the unit normal vector n of the triangle X{A), 
and (b) the unit normal vectors iij (i = l,2,3) of the nearest neighbor triangles. 



The partition function of the model on fixed connectivity surfaces, which is 
denoted by model 1, is defined by 



ZUb) = I [dL] I [dip] [ n dX, exp [-5(X, {L, if})] , (model 1) 

i=l 



(9) 



where the symbols J[dL] and J[d(f] are given by 



and 



'[dL] = ll[dL,expi-J2LU: 

i=i \ 1=1 J 



(10) 



The symbol Li in Eq. ffTOj) denotes the bond length of A, and ipi in Eq. ffTTl) 
denotes the deficit angle. The factors exp{—J2iLf) and exp {— are 
necessary in order to make the integrations of the variables L (g R>o) and 
(/^(gR) well-defined; the factor exp {—J2a ^a) can also be used in place of 
exp{—J2iLj) in Eq. ffTUl) . It is possible to introduce parameters A^, in 
order to control the fluctuations of L and (p such that exp (^—Xl J2^=i -^i) and 
exp (^—X^Yli=i however, we assume = 1 and A^ = 1 for simplicity. 

The prime in /' Hili dXi in Eq. denotes that the center of mass of surface 
X{M) is fixed. We note that M is not always globally fiat (or MCR^) even 
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when $ = $0 because the deficit angles at the vertices are not always zero, this 
is because the edge length L is integrated independently of the other edges. 

The partition function of the model on dynamically triangulated surfaces, 
which is denoted by model 2, can be written by including the sum over possible 
triangulations J2t such that 

Z{iuib)=Y. f[dL] [[dip] [ f[dX,exp[-S{X,{L,ip})], (model 2). (12) 

T i=l 

X^T is performed by the bond flip technique. The bond flips is simultaneously 
performed both in M and in X{M), because the mapping X preserves the 
triangulations. Since the edges of triangles in M are considered as coordinate 
axes, the flip of bonds can make a large difference on the configuration and 
influence the equilibrium property. Thus, the bond flips should be carefully 
performed in model 2. The detailed information on this point is given in Sec- 
tion H 



3 Continuous model 



We comment on a correspondence between the discrete model and the con- 
tinuous model. The continuous energy Si is the Polyakov action for strings in 
R'^ {d = 2), and it is given by 

^1 = / Vad^xg'^'daX'^dkXf^, (13) 

where g is the determinant of the metric tensor gab, (a, 6= 1, 2) of M, and g"''^ 
is the inverse of gab- The symbol fi of X^ denotes that X^ G IV^ {d = 2). In 
order to obtain an explicit expression of gab, we consider the edges Li and 
L2 of Z\ in M as the axes of a local coordinate. Thus, we define the discrete 
metric gab such that 



Sab 



^ L\ LiL2Cos$3 
^LiL2Cos$3 L\ 



(14) 



We should note that gab is not exactly identical with the induced metric of the 
co-ordinate mapping \& from A to R"^ {D = 2) but gab is close to the induced 
metric of \E'. In fact, gab is just identical with the induced metric of \1/ if 
$ = $0- We note also that gab is independent of X, and it depends only on 
L and ip; both L and p) are considered as functions on A. If gab is given by 
the induced metric of X such that gab = daX^dbX^^, then A (c M) can be 
considered as a subspace of IV^{d = 2), and L and $ are identified with i 



7 



and respectively, and the area is just identical to the area of X{A). In 
this case, the length scale of M is automatically determined by the Euclidean 
length scale of R*^ {d = 2), and hence the length scale of M remains unchanged 
no matter how X{A) and as a consequence g transforms its shape according 
to the variation of X. We note also that gab in Eq- (|T4l) makes the form 
ds^ = gabdxadxh positive definite. 

The partial derivatives daX'^ are replaced by X2 — X1 for a = 1 and X3 — X1 
for a = 2, where Xi and X2 are two terminal vertices of the bond ii in the tri- 
angle shown in Fig. [U^c), and Xi and X3 are those corresponding to the bond 
£2- Since / ^/gd'^x corresponds to the area of M, / y^rf^x can be replaced by 

Aa- Thus, we obtain ^1 = ^4 (1/4A^) {Llil+Llil-2iii2 cos 03L1L2 cos $3) 
Including the terms, which are the cyclic permutations such that 1— )-2, 2— >-3, 
and 3 — )■ 1, and using the multiplicative factor 1/3, we have the expression 
of Si in Eqs. (|6]) and ([7]). This symmetrization makes 5*1 reparametrization 
invariant in the sense that three pairs of edges (Li, L2), {L2, L^), and (L3, Li) 
can be considered as the coordinate axes. 

If the variable L and $ in Si of Eq. (jH]) are replaced by i and respectively, 
then Si{i, (p) coincides with twice the area of X[M) as a subspace of R'^ ((i = 2); 

0) =2 X]z\ Q-zi, where aA is the area oiX{A). If in Eq. ( !T3|) is multiplied 
by the factor 1/2 such that (l/2)S'i, then the factor 1/12 of 5*1 in Eq. ([6]) is 
replaced by 1/24, and as a consequence the corresponding discrete 5*1 (£,0) 
is just identical to the area of X{M). We should note also that the phase 
structure of the model is independent of the multiplicative factor such as 1/2 
of ^1. 

The continuous bending energy 5*2 is given by 

S2 = \j ^d^xg'^'dau^dbu^, (15) 

where is a unit normal vector of the continuous surface in R*^ {d = 2) and 
has values in Z2 = {1, — 1}; the symbol /i of can be dropped. We should 
note that the expression of 5*2 of Eq. ( |T5l) coincides with that of Polyakov's 
extrinsic curvature term (1/2) / ^/gd'^xK^K^ in if c? = 3 and gat is assumed 
to be the induced metric such that gab = daX^dbX^, where Kab is the second 
fundamental form defined by Kab = —daX^^dbU^ ■ In fact, it is straightforward 
to see this by using the equality dati^ = —K^dbX^^. We should note that S2 
of Eq. fll5p is not always identical to the Polyakov's extrinsic curvature term, 
because gab in Eq. ffT^ is different from the induced metric of X. 

The normal vector n'^ is defined on the triangles, and therefore we replace din'^ 
and d2n'^ such that c?in^— )-no — n2 and c?2n'^— J-nQ — rii, where iij are shown in 
Fig.[2](b). The discrete version of S2 is then given by X]zi(l/4^/i)[-Z^i(l— no-ni)+ 
L|(l— no-n2)— I/1L2 cos $3(111— no){n2— no)]. Symmetrizing this term by the cyclic 
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permutations just like in the case of 5*1, we have 5*2 = (1/12) X^zi 'S'2(/i)/^/i, 
where 5*2 (Zi) is given by Eq. ([8]). 

We should comment on why 5*2 in Eq. f|T5|) is multiplied by the factor 1/2, 
which makes the factor of the discrete 5*2 in Eq. ([8]) as 1/12. This is because 
5*2 in Eq. is convenient to compare the results with those obtained from 
the conventional model, whose definition will be mentioned below in Section 
[51 Because of the factor 1/12, the value of S2/NB is comparable to the one of 
the conventional model at the same b. 

The continuous partition function is expressed by 



where / DX denotes the sum over mappings from M to R'^, and / Dg the 
sum over metrics on M. The integrations / DX and / Dg are considered to 
represent the sum over surfaces in the string model context. However, we are 
not going into details of the measures [24j. In this paper, / DX and / Dg are 
simply replaced by the three-dimensional multiple integrations /' Hi^i dXi in 
Eq. IQ and J[dL] J [dip] in Eqs. (dO]) and (III]), respectively. 



4 Monte Carlo technique 



The vertex position X (^E R*^ {d = 2)j is moved to a new position X' = X+6X, 
where SX is chosen randomly in a small sphere. The radius of the small sphere 
is fixed as an input parameter of the simulations. The edge length L and the 
deficit angle if are also changed to L' = IA-5L with L' > and ip' = (f-\-6(p, where 
6L and 6(p are chosen randomly in small one-dimensional ranges [—Lq, Lq] and 
[—ipo, (po]. The constants Lq and (po are fixed to Lo = 0.5 and ipo = 0.25. Triangle 
equalities of Eq. ([T]) are also imposed on the variation of L. 

Both of the energies Si and S2 vary not only with 6X but also with 6L and S(p, 
because the area of the triangle A depends on edge length L and the inter- 
nal angle $ at least. The internal angle $ of Z\ is also dependent on the deficit 
angle ip. Three variables cp, L, and X are assumed to be independently varied; 
a variation of one variable does not change the remaining variables. Moreover, 
the intrinsic variables ip and L of one triangle are assumed to be independent 
of those of the other triangles in M. Only constraint on the triangles is the 
fact that L is a common variable shared by two neighboring triangles. These 
assumptions are understood as reasonable because the variables can be consid- 
ered as functions on Z\ C M. We should note that the assumptions are illegal 
if g is given by the induced metric gab = daX'^dbX^, where M is considered to 
be included in R*^ {d = 2), and in this case the bond length of linear triangles 




(16) 
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varies when the vertex position varies, and as a consequence the length scale 
of M is fixed to the Euclidean length scale of the bulk space R'^((i = 2) as 
mentioned in the previous section. 

The triangulated surfaces are obtained from the surfaces constructed in R'^ ((i = 3) 
by assuming the third- component of the vertex position as zero. The surfaces 
in R'^((i = 3) are obtained from the icosahedron by splitting the edges into 
small pieces and are identical to those used in [21] . 

One Monte Carlo sweep (MCS) consists of updates of X, A''^ updates 
of L, and Nt updates of ip. Not only the vertices but also the bonds and 
the triangles are labeled by sequential numbers, and therefore the Metropolis 
updates are done by using these sequential numbers. The new values of the 
variables are accepted with the probability Min [1, exp(— 55")], where S' is the 
effective Hamiltonian including the terms from the integration measures in 
Eqs. (IIOD and ^ such that S' = Si+hS2+Y.i=iL'i+T.7=i We have about 
50% acceptance rate of X, 78% acceptance rate of L, and 93% acceptance rate 
of Lp. The acceptance rates of L and are almost independent of variations of 
Lq and ifQ. 

The bond fiip technique is assumed to define the sum over triangulations X^t 
in the case of fiuid surfaces. Flip of bonds is performed both on M and X{M) 
in 'Rf{d = 2). The coordination number q is bounded such that 3 < g < 30. 
The phase structure seems not to be so strongly infiuenced by the assumed 
upper bound gmax = 30, because almost all q are smaller than gmax = 30 in 
the configurations at relatively small h. The length of fiipped bond on the 
surface X{M) is automatically determined by using the canonical coordinate 
of R"' ((i = 2), while the length L' of the fiipped bond in M is randomly chosen 
such that L' = L + 6L, where L is the length of the bond to be fiipped and 
6L{e [—0.5,0.5]) is a random number. The reason why L' is given by such 
definition is because we have no information on the length of fiipped bond in 
M. Information on M is obtained only locally and limited to a triangle A and 
its tangential triangle; the bond length L shared by two triangles only connects 
one triangle to the other on M. This definition of the fiipped bond length 
seems to be a reason why the bond fiip strongly infiuences the equilibrium 
configurations as mentioned in Section [21 Thus, we perform the bond fiip 
N/npinp = 100) times a MCS by choosing bonds randomly. If the fiips are 
performed more frequently; for example N/npinF = 2) times a MCS, the 
expected relation Si/N=l is considerably violated. This relation is expected 
due to the scale invariance of the partition function, and no violation is seen in 
the model defined without the variables L and even when the fiip of bonds 
is performed more frequently. The rate of acceptance of the bond fiip is about 
40% ~ 50%, which is almost independent of the magnitude of 6L, and the 
rates of acceptance of the remaining variables are almost the same as in the 
case of the fixed-connectivity model. 
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In order to see the dependence of the results on N/np, we perform the simula- 
tions on the conventional fluid model under N/np = 5, which is different from 
the condition N/np{np = 100). The condition N/nplnp = 100) is given by fix- 
ing Up to be independent of the size N, while the condition N/np = 5 is given 
by varying np to be dependent on such that np = N/5. We find that the 
critical exponents are almost independent of the conditions. The simulations 
with N/np = 5 is relatively time consuming than those with N/np{np = 100), 
because the convergence speed is relatively low, or in other words a large num- 
ber of MCS is necessary to obtain high statistics data, in the simulations with 
N/np = 5. This is because the total number of bond flips (=5) in one MCS is 
very small in this case. Thus, the condition N/np{np = 100) is assumed also 
in model 2. 

The total number of MCS for the production runs after sufficiently large num- 
ber of MCS for the thermalization is 6 x 10^ ~ 9 x 10® at the transition re- 
gion of the = 2562, = 3612 surfaces for model 1, and relatively small 
number of MCS is assumed at the non-transition region and on the smaller 
surfaces. Almost the same total number of MCS is assumed in the conven- 
tional fixed-connectivity model. For model 2 and the conventional fluid model, 
6 X 10® ~ 9 X 10® MCS are assumed at the transition region of the A^ = 2562, 
A^ = 3612 surfaces, and relatively small number of MCS is assumed at the 
non-transition region and on the smaller surfaces. 



5 Numerical results 

5.1 Fixed- connectivity surfaces 

In this subsection, we show results of the fixed-connectivity model, which is 
model 1, and those of the conventional model defined only by the variable 
X (^G R'^((i = 2)^. The partition function of the conventional model is given 
by 



^c.fix(&)= / l[dX^exp[-S{X)], S = Si + bS2, (17) 

i=l 

Si = Y. (^^ - ' ^2 = E (1 - ■ n,) , n, G Z2 = {1, -1}, 

where Si and 5*2 are the Gaussian bond potential and the bending energy, 
which correspond to Si and 5*2 in Eq. ([6]). The symbol c.fix in Zc.fix(&) denotes 
the conventional fixed-connectivity model. 
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b 



Fig. 3. (a) The mean square size X'^ vs. b of model 1, (b) the variance C^a vs. b, 
and (c) the peak values C^2^ vs. in a log-log scale. The solid curves in (a) and 
(b) are drawn by the multi- histogram reweighting technique |28j . 

Figure W^a) shows the mean square size X"^ vs. b, where is defined by 

1 Af _ 9 _ 1 ^ 

^^ = ^E(^^-^) , ^ = ^E^^, (18) 

where X is the center of mass of the surface. The sohd curves in Figs. [3]^a) and 
^h) are drawn by the muhi- histogram reweighting technique [28]. We see that 
X"^ smoothly varies against b, and that rapidly increases with increasing 
N. On the surfaces of N <2562, X"^ changes up and down many times during 
the simulations at the transition point. To the contrary, on the largest surface 
of A^ = 3612, the surface configuration seems to be trapped in one of the two 
potential minima at the transition point. Thus, the transition is not always 
correctly refiected on the X = 3612 surface. 

The variance Cx^2 of defined by 

Cx'^ = l^{{x'-{X')y) (19) 

is plotted in Fig. Mjo)- We find that Cx^ has a peak, and the peak value C^f""" 
increases with increasing N. This indicates that the shape transformation is 
refiected in the fiuctuation of X^. 

The peak values C™2^ are plotted in Fig. [31(c) against X in a log-log scale. The 
straight line is drawn by fitting the data to C™2^ ~ X°", where C™2^ of the 
X = 2562 surface is the value obtained by the multi-histogram reweighting. 
Thus, we have 

cTi = 1.77 ± 0.07, (model 1) . (20) 

We examined the exponential fitting such that C™f^ ~ exp{a), however, the 
power law fitting C^?^ ~ X"^ is better than the exponential fitting. This 
observation can also be seen in all other physical quantities in all of the models, 
which will be studied in this paper. 
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The result ai = 1.77(7) in Eq. ( 120|) clearly indicates that the transition is of 
first order. 




Fig. 4. (a) The mean square size X'^ vs. 6, (b) the variance C^a vs. b, and (c) 
the peak values vs. A'' in a log-log scale. The solid curves in (a) and (b) are 

drawn by the multi-histogram reweighting technique. The data are obtained by the 
conventional fixed-connectivity model of Eq. (jlTp . 



Figures lll^a)-[4l^c) show the results of the conventional model of fixed connec- 
tivity surfaces defined by Eq. ( flTl) . The data shown in the figures correspond 
to those shown in Figs. [3](a)-[3](c). The phase transition is seen in the iV = 3612 
surface of the conventional model in contrast to the case of model 1. In fact, 
the surface configuration seems not to be trapped in one of the potential min- 
imum states at the transition point; this can be seen in the variation of X 
against b, and for this reason the variance Cx^ is correctly computed at the 
transition region even on the = 3612 surface. To the contrary, as we see 
in Figs. El^a) and [3](b) the surface configuration seems to be trapped in the 
potential minimum states in model 1 on the A^ = 3612 surface. 



So/N 



(b) I model 
N=3612 

N=2562 




0.49 0.5 , 0.51 

b 



0.49 0.5 , 0.51 

b 




1000 2000 3000 



Fig. 5. (a) The bending energy S2/NB vs. b, (b) the specific heat vs. 6, and (c) 
the peak values Cg]^^ vs. in a log-log scale. The data are obtained by model 1 of 
Eqs. dSD-dllJ. 



The bending energy S2/NB and the specific heat are plotted in Figs. [5] 
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Fig. 6. (a) The bending energy S2/NB vs. b, (b) the specific heat vs. b, and 
(c) the peak values Cg]^^ vs. in a log-log scale. The fitting was done by using the 
largest four data in (c) . The data are obtained by the conventional fixed-connectivity 
model. 

and El where is defined by 

Cs2=^^{{S2-{S2)f). (21) 

The results of model 1 are shown in Fig. [5], and those of the conventional 
model are shown in Figs. [61 The data C^*" of the N = 2562 surface in Fig. 
[5]^c) is the result of the multi-histogram reweighting, since C^^**^ is slightly 
smaller than the peak of the solid curve as we see in Fig. [51(b). The straight 
line in Fig. M^c) is the fitted one of data to Cg"^"^ ~ A^": 

tti = 1.98 ± 0.07, (model 1) . (22) 

The data C^^^ and the fitted line of the conventional model are shown in 
Figs. [61(c). We see that S2/NB of model 1 appears to be trapped in one of 
the two different values on the = 3612 surface at the transition point in 
contrast to the conventional model, where S2/NB changes up and down many 
times during the simulations on the A^ = 3612 surface at the transition point. 
However, the phase structure of model 1 is considered to be identical with 
that of the conventional model, because the value of ai of model 1 is almost 
exactly identical with a of the conventional model. 

Finally in this subsection, we plot the Gaussian bond potential Si/N vs. h in 
Figs. [71(a) and[71^b). From the scale invariance of the partition function, Si/N 
is expected to be Si/N= {N — 1)/A^~1. We see from the results in Figs. [Tl^a) 
and [71(b) that this expectation is satisfied in both models. 
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Fig. 7. The Gaussian bond potential Si/N vs. b of (a) model 1 of Eqs. ([6|)- (fTT]) and 
(b) the conventional fixed connectivity model of Eq. (jl7p . 



5.2 Fluid surfaces 



The fluid surface model denoted by model 2 is defined by the partition func- 
tion of Eq. f|T2|) . which includes the sum over triangulations. Flips of bond 
discontinuously change the surface configuration in contrast to the cases of 
continuous variations of X, L and if, and hence, the bond flip is performed 
only N/np{nF = 100) times a MCS as described in Section HJ The results are 
compared with those of the conventional fluid model defined by the partition 
function 

/ N 

Zc.&n{b) = E / n exp [-S{X)] , S = S, + bS2, (23) 

T i=l 

where Zc.an denotes the partition function of the conventional fluid surface 
model, and 5*1 and 5*2 are given by Eq. ( fTTj) . The total number of bond flips 
N/np per one MCS is also assumed to be N/nplnp = 100) in the conventional 
fluid model. The length of flipped bond is automatically obtained in the con- 
ventional model, and hence it seems not necessary to reduce N/np so small, 
because the equilibrium surface configuration is not so strongly influenced by 
the bond flip in contrast to the case of model 2. However, the same np is 
assumed in the conventional model as that of model 2 in order to compare the 
results under the same condition. 

Figures [8] and |9] show the mean square size and the variance C^z of model 
2 and the conventional model. The peak values C^f^ are shown in a log-log 
scale. We find that the transition is of second-order in both models. The fitted 
value of the critical exponent, which is defined by C™2^ ~ A^'^, is 

£72 = 0.86 ± 0.04, (model 2) . (24) 

The value of cr2 is (72 < 1, and this implies that the transition is of second 
order, although it is close to first order because cr2 — 1. Thus, the order of the 
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Fig. 8. (a) The mean square size X"^ vs. b, (b) the variance C^a vs. b, and (c) the 
peak values C^2^ vs. in a log-log scale. The solid curves in (a) and (b) are drawn 
by the multi-histogram reweighting technique. The data are obtained by model 2. 
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Fig. 9. (a) The mean square size vs. b, (b) the variance C^a vs. b, and (c) the 
peak values C™2^ vs. in a log-log scale. The data are obtained by the conventional 
fluid model of Eq. (p3|) . The data denoted by symbols Q ^ ™ (c) correspond to 
the simulation conditions iV/rij?(ni? = 100) and N/np = 5, respectively. The fitting 
was done by using the largest four data in (c). 

transition remains unchanged, although the exponent a2 of model 2 is slightly 
larger than cXc.flu. = 0.74(6) of the conventional model. The intrinsic variables 
L, (fi slightly strengthen the transition of fluid surface model in contrast to the 
case of the fixed connectivity model in the previous subsection. In Fig.[9]^c), we 
show the data denoted by the symbol (A), which are obtained with the simu- 
lations under the condition N/nF = 5. The exponent CTcfiu. = 0.76(6) is almost 
identical to the cXc.flu. =0.74(6) obtained under the condition N/nF{nF = 100), 
although Cx2^ slightly depends on the conditions. This implies that the final 
results are independent of the simulation condition of N/np in the conven- 
tional fluid model. Thus, it can also be expected that model 2 is independent 
of the condition, because N /npirip = 100) is considered to be sufficiently small. 



The bending energy S2/NB and the specific heat of model 2 and the 
conventional fluid model are shown in Figs. [10] and [TTl where is defined 
by Eq. (12T!) . The peak values C"^^^ grow larger with increasing iV, and this 
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Fig. 10. (a) The bending energy S2/NB vs. b, (b) the specific heat vs. b, and 
(c) the peak values Cg^^^ vs. in a log- log scale. The fitting was done by using the 
largest four data in (c). 
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Fig. 11. (a) The bending energy S2/NB vs. 6, (b) the specific heat vs. b, and (c) 
the peak values C^^^ vs. in a log-log scale. The data denoted by symbols Q and 
A in (c) correspond to the simulation conditions N/npiriF = 100) and N/np = 5, 
respectively. 

behavior of C^*" of model 2 shown in Fig. [TOT b) is almost identical to that of 
the conventional model in Fig. [TTT b). Cg^^^ is plotted in a log- log scale against 
N in Figs. [TUT c) and [TTT c). The exponent a defined by C^^^ ~ A^", which is 
the slope of the log- log fit, is given by 



02 = 0.72 ± 0.07, (model 2) 



(25) 



where the fitting was done by using the largest four data. The value of 02 is 
slightly larger than a = 0.65(5) of the conventional model. This observation is 
also consistent with the previous ones that the order of the phase transitions 
of model 1 and model 2 remains unchanged from those of the conventional 
models. In Fig. [TlT c). we show the data (A) obtained with the simulations 
under the condition N/np = 5. We find that the exponent a = 0.63(5) is almost 
identical to a = 0.65(5) obtained under the condition N/np{np = 100). 



Figures [T2](a) and [TW b) show the bond potential Si/N vs. b. The potential 
Si/N in Fig. [TWb) is slightly lower than Si/N = 1, this is because Si/N = 



17 



Si/N 
1.002 

1 

0.998 



(a) model 2 




^7 ^7 



8= 



J I I 1 I 1 L. 



Si/N 



1.002 



0.998 




a;N= 


3612 


o:N= 


2562 


0:N= 


1962 


□ :N= 


1442 


v:N= 


1002 



1 1.02 1.04 , 1.06 

b 



0.91 0.93 , 0.95 

b 



Fig. 12. The bond potential Si/N vs. b of (a) model 2 and (b) the conventional fluid 
model. Si/N has the expected value S'i/A^ = l. 

(A^— 1)/A^ just as the one in Fig. [71 In the conventional model, the length of the 
flipped bond is exactly obtained, because the triangulated surfaces is included 
in R^. Therefore, one can expect that the equilibrium configurations are not so 
strongly violated by the bond flips. To the contrary, the length of flipped bond 
in the triangulated surface in M is randomly chosen as described in Section 
m Thus, bond flips can influence the equilibrium property of configurations 
of model 2. This is a reason why Si/N in Fig. [T^ a) is slightly larger than 
Si/N = 1. In fact, the deviation of Si/N from Si/N = 1 grows larger when 
N/np the total number of bond flip per one MCS is assumed to be N/np > 
A^/100. If Up is assumed to be np > 100, then we have Si/N which is more 
close to Si/N=l. The deviation of Si/N seems to grow with increasing in 
FigJTWa). This implies that np should be increased with increasing N. 



Finally, we show the deficit angle (p vs. b in Figs. [T3l(a) and fTST b). The symbol 
(fi denotes J^i^Ji^fi/^T, where (pi is defined by Eq. ([3]). The variations of (p 
against b are similar to those of S2/NB in both model 1 and model 2. The 
discontinuity seen in (p of model 1 is very small compared to the value of ip 
itself, however, we see that the phase transitions are clearly reflected in the 
internal geometric variables. The variance defined by 



Nt 
i=l 



(26) 



is plotted in Figs. \T3\ c) and [T3l^ d). The shape of is similar to those of Cx^ 
and in each model, however, the peak value C^'^^ increases only slightly 
with increasing N. For this reason, the scaling of the peak value C™^^ such 
as C^'^ ~ A^^ is observed neither in model 1 nor in model 2. Although the 
phase transitions of the models are reflected in ip, the gap of ip is not always 
considered as a signal of a transition in M. 

We should note that J2i=i P^i/Np is expected to be zero in the limit of A^— J-O 
on smoothly triangulated surfaces. Thus, non-zero p aX N ^00 is possible in 
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Fig. 13. The deficit angle i-p vs. h of (a) model 1 and (b) model 2, where ip is defined 
by 99 = YL^Zx "Pi I^T, Nt = 2N-4. The unit of if is [degree] , which is [radian x ISO/vr] . 
The variance is shown in (c) and (d), where no error-bar is shown in (d). The 
values of C^{ip = Nt y:>) in (c) and (d) are reduced by a factor (vr/lSO)^ if the unit 
of if is changed from [degree] to [radian] . 



the model of this paper as mentioned in Section [21 and in fact it is clear that 
(p>0 in both model 1 and model 2 at the transition points at least. 



We finally comment on the local fluctuation of (fi, which is not presented 
as a figure. In the simulations on both fixed-connectivity and fluid surfaces, 
the minimum v'f™ and the maximum (yj™^^ are respectively comparable to 
±27r(~ iSvr), which are considered to be out of the range — 7r<(/?<27r, which 
corresponds to < $ < vr. Moreover, the mean value of \ip\ is about 7r/3.5 in 
both fixed-connectivity and fluid models, and therefore the local fluctuation 
of (fii is very large compared to the mean value of if estimated from the data in 
Figs. [TST a) and lTST b). The large local fluctuations of ipi seems due to the fact 
that no interaction of (fi is assumed in the model of this paper. However, we 
consider that there is no influence of such a relatively large local-fluctuation 
of if at least on the phase structure as we have confirmed from the presented 
numerical data. 
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6 Summciry and Conclusion 



We have numerically studied a triangulated surface model, which is defined by 
a mapping X from a two-dimensional spherical surface M to 'Rf' {d = 2). The 
dynamical variables of the model are the metric g oi M and the mapping X, 
which are summed over in the partition function. Hamiltonian S of the model 
is given by a linear combination of the Polyakov action Si for strings and the 
extrinsic curvature S2 such that S = Si + bS2, where b is the bending rigidity. 

By using the Regge calculus technique, the integration over g in the partition 
function is replaced by the integrations of the edge length L and the deficit 
angle ip of the triangle A in M. The variable g is defined to be a small variation 
of the induced metric of the coordinate mapping from AinM toR^ {D = 2), 
where the variation is given by the deficit angle ip. If ip is assumed to be 93 = 0, 
g is just identical to the induced metric of ^. In this case, M is still not 
always completely flat even though A becomes a hnear triangle. Thus, S is 
defined to be dependent not only on the extrinsic variable X but also on the 
intrinsic variables L and ip; S = S{X, {L, ip}). The integrations of the variable 
X and the variables L and ip arc performed in MC simulations by deforming 
the triangulated surface X{M) in R'^ {d = 2) and the triangulated surface M, 
respectively. We should note also that the triangle inequalities are strictly 
satisfied not only on the triangles in X{M) but also on the ones in M during 
the MC simulations. 

Our attentions are focused on whether the intrinsic variables influence the 
phase transitions corresponding to the surface fluctuations and the collapse 
phenomenon. In order to see this influence we study the two variations of the 
model; the flrst is the flxed-connectivity model, and the second is the fluid 
surface model, which is deflned on dynamically triangulated lattices. Since 
the triangle edges on M play a role of local coordinate axes, the dynamical 
triangulation is considered to make the model reparametrization invariant. 
The conventional model, which is deflned only by using the variable X, is also 
studied in order to compare the results with those of the models in this paper. 

Our conclusion is that the internal geometry does not so strongly influence 
the transition of shape transformation. The order of the transition is of flrst 
order in the flxed-connectivity model and of second order in the fluid model. 
The order of the transition remains unchanged from the corresponding con- 
ventional model on both flxed-connectivity and fluid surfaces, although the 
critical exponents of the transition are slightly different from each other in the 
case of fluid model. It is also found that the deflcit angle (p (^Yji'=i V^i/^t) dis- 
continuously changes at the transition at least in the flxed-connectivity model, 
and hence the transition is reflected in internal geometric variables. 
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It is interesting to study the surfaces embedded in R*^ {d = 3). It is also inter- 
esting to study the case where Qab depends only on the variable (f, and the 
case where gab is not always given by an induced metric such as the one in 
this paper. These remain to be studied in the future. 
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